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Abstract 

A simple formula for correlation energy E c of the ir electron systems is obtained under an approximation for the electron- 
electron interactions. This formula is related directly to square of the bond order matrix and the nearest-neighbor Coulomb 
electron-electron interaction. The influence of the correlation energy on the band energy gap is discussed. The values of the 
correlation energy for polyacetylene (PA) are calculated and can be compared with those for PA obtained by other methods, 
including ab initio method. 
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£ Introduction 

o : 

i ^rne electron correlations have been a very important issue in investigating the electronic structures of various electron systems. 
Especially the electron correlations have a strong influence on the bonding properties of atoms and semiconductor band gap 1 . 
CftWever, it is well known that it is very hard to completely solve a many electron system analytically using a single approximation 
^jcause of the exchange and correlation problems. Hartree-Fock approximation deals with the exchange problem between the 
j^eune spins among electrons but not resolve the correlation between the opposite spins among electrons 2,3 . The local density 
^approximation (LDA) of Kohn and Sham 4 and later the density functional theory of Hohenberg and Kohn 5 made a contribution 
l/o) the exchange-correlation energy, denoted by E ex [p(r)\ using a complicated functional form. In the LDA, the ground state 
^sfctiange energy E x can be expressed an integral of the charge density. However, it is difficult to express the correlation energy 
jSr: in an available form of integral of the charge density which is easy to calculate and only could be expressed as a numerical 
■~S3rmula with the parameter r s after a large number of works 6 . 

Besides the LDA, there are other methods which can be used to deal with the correlation effects of various electron systems 
Sjjch as metal and nonmetal atoms, small and big molecules, as well as polymers. They include unrestricted self-consistency 
hjeld (SCF) Approximation 3 , Configuration Interaction method(CI) 7 , coupled-cluster method 8,9 that are applied to quantum 
' fcJatemistry and nuclear physics 10 , Jastrow wave-function method 11 that is used to describe correlations in homogeneous fermion 
Stems like the electron gas or liquid He-III 12,13 , the Projection Technique 14 for the strongly correlated systems, and the GW- 
approximation (the one-particle Green's function plus the screened Coulomb interactions) by Hedin 15 . For the conjugated 
polymers such as polyacetylene (CH) X (PA) or poly(p-phenylenevinylene) (PPV), the exchange-correlation energies has been 
.Peculated by the method of the ab initio with GW-approximation 16 . For the various metal and nonmetal atoms, one used 
^s^omic Bethe-Goldstone equation under Hartree-Fock functions to calculate the electron correlation energies for their ground 
Hates 17 > 18 < 19 . For small or medium-size molecules(hydrocarbon molecules), a semi-empirical SCF scheme plus CNDO or INDO 
approximations was used to calculate inter-atomic correlations and intra-atomic correlations whose results could be compared 
with ab initio method 20 . 

However, usually the discussion about the correlation energy concentrates mostly on the on-site Coulomb interaction, that 
is, from the term m\nii and seldom specially on the lang-range Coulomb interaction, at least, on the nearest-neighbor Coulomb 
interaction. For strongly correlated systems, the on-site Hubbard interaction U is much bigger than the nearest-neighbor 
Coulomb interaction v, so the contribution of the correlation energy are mainly from U . But it is well known that in the most 
conjugated polymers, the on-site electron-electron interaction U is not so bigger than the nearest-neighbor interaction v because 
of screening 21 . 

In the past, people treated with the correlation energies both from the on-site Hubbard interaction U and the long-range 
Coulomb interaction, In Ref.22, authors studied the correlation energies of polyethylene (CH2) X (PE) using the local ansatz 22 . In 
Ref.23, authors used the Gutzwiller ansatz as the variational ground state and studied correlation energy of polyacetylene (PA). 
First they used the Hubbard term plus SSH Hamiltonian to discuss the correlation energy and later added the nearest-neighbor 
interaction to get an effective U e ff and discuss the correlation energy again. In Ref.24, author used a variational method 25 ' 26 to 
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study the correlation of PA in the PPP model where both Hubbard term and the long-range Coulomb interaction were included. 
But author did not specifically deal with the correlation energy from the long-range Coulomb interaction. It is important for 
the long-range interaction to calculate the band energy gap in the conjugated polymers 27 . 

Although there were those studies above about the correlation energy both from the on-site Hubbard interaction U and 
the long-range Coulomb interaction and people have been known the importance of the correlation between two tt electrons in 
understanding electronic properties of the conjugated polymers, such as the optical band gap, exciton excitation energy and 
exciton binding energy, polaron, etc., there is still lack of special studies on the correlation energy mainly from the long-range 
Coulomb interaction. Therefore, study of the correlation effect due to the Coulomb interaction of two adjacent it electrons in 
the conjugated polymers becomes significant. 

The purpose of the present work is to study the correlation energy mainly from the long-range Coulomb interaction (the 
nearest-neighbor interaction in this study) . The starting point in this paper is the exchange-correlation energy E ex . Although the 
correlation energy may be expressed as E c — J e c n(r)dr, e c is hard to know and hard to obtain. Thus, available approximations 
will be adopted in this study. Finally, a formula of the electron correlation energy due to the nearest-neighbor Coulomb 
interaction is obtained in an analytical form. Using this formula, it is simple to calculated the correlation energy of a ir 
conjugated polymer and available to discuss the influence of the correlation energy on the band energy gap. 

The arrangement of this paper is as follows. In section II, the expression of the pair-distribution function for the conjugated 
polymer systems by omitting wave function overlap integrals between two adjacent atomic sites is presented. In section III, an 
approximation for the electron-electron interaction integral is made and an analytical expression of the correlation energy for 
the conjugated polymers is obtained. In section IV, the correlation energy calculation for one-dimensional polyacetylene (PA) 
chain under tight-bind approximation (SSH Hamiltonian plus electron-electron interactions), and results are presented, and in 
section V is discussion with a summary. 



2 Pair-distribution function for electron systems 

The so-called correlation is the correlation hole appearing around an electron moving in the medium. Remarkably the correlation 
in fact is due to the electron-electron interaction among electrons. In a homogeneous or an inhomogeneous electron gas, the 
density operator of the electrons is expressed as a delta function p{r) — ^2n$(r — r*j), the density distribution function of the 
electrons is the average of the operator over a given state |$ >, p(r) =< Q\p(r)\$ >. The two-point density-density correlation 
function at a given state |$ > is defined as 

C(f,P) = < <l>\^26(? -ri)5(f-fj)\<f> > 

= <p(r)p(f>)> (1) 

where < . . . > means < $ . . . |$ >, the average of the operators over the ground states. To calculate the two-point density- 
density correlation function, one introduces the pair-distribution function which is defined as 28 

< p(f)p(?) >= g(r,r') < p(r) >< p(?) > (2) 

About the pair-distribution function g(f,r') there are different approximations. For example it can be cast into an analytic 
expression for the homogeneous gas 29 . For a many-electron system, the density distribution p(f) of the electrons can be expressed 
by 

occ 

p(r) = 2^|^(r)| 2 

N 

= EptftfWiW (3) 



where the index p refers to the occupied molecular orbital, p = 1, 2, • • • , N/2. The indices i and j refer to the atomic sites. 4>i{r) 
is the electron orbital wave function at the atomic site i and is normalized to one. Here pij = < a/ ia a,j > where a\ a (a,j a ) 
is the electron creative (annihilation) operator with spin a at the site i (j). pij is called the bond-order matrix 28 . Comparing 
(3) with p(r) =< <j>|/5(r)|$ >, we can see that the density operator for a many-electron system can be represented by 



N 



Kr)=J2P^*^(r) (4) 
with pij = J2 a a la a jv Then the two-point density-density correlation function is expressed as 

N N 

< p(r)p(P) > = EEE < P°jCi>UVMVril<t>Ur')Mr') > 

aa' ij kl 



N 

= E E <PuPik>\Mr)\ 2 \Mr')\ 2 

<x<x' ik(i^k) 

N N 

+ EE E <PijPia>4>mMWk&)MJ) (5) 

and 

N 

<p(r)> - <E^«(%«> 

TV 

= E<^>i^i 2 +EE<^ > «('i (6) 

A calculation containing the pair-distribution function involves wave function integrals between different sites, which is related 
to the overlap of two electronic wave functions. As a first order approximation and without loss generality, the overlap are not 
considered, so the second terms in (5) and (6) could be omitted in this study. Therefore, the pair-distribution function g(f,r') 
in (2) becomes 

? = Ef fc( ^ ) <^;>i^^i 2 i^(^)i 2 

E CTi < Pit > MWZrk < Plk > IMr'W 

with g(f,r') = J2*v> 9(r,r') and p^ =< a\ a a ja >. 

3 Approximation approach to the correlation energy 

In the local density approximation (LDA) the exchange-correlation energy can be written, in terms of the pair-distribution 
function g(f,r') in the form 30,31 



where 



E XC [P\ = ^ I dW P (r) W l r,) l] p{?) 
^ J \r — r\ 

g(f,r') = / d\g(f,P;X). 
Jo 



(8) 



(9) 



Here g(r,r') includes the exchange-correlation information between two spatial points r and r'. Here A means the charge e 2 in 
Coulomb interaction is replaced by Ae 2 in the process of calculation. In (8), p{r')[g{r, r') — 1] describes the depletion in density 
due to the exchange-correlation hole around an electron at position r. The density depletion due to the exchange-correlation 
hole around an electron corresponds to exactly one particle, independent of the size of the electron-electron coupling strength. 
Then we have the sum rule satisfied by the exchange-correlation hole: 



J d 3 r'p{r')[g(r,r') - 1] = -1. 



(10) 



It is well known that it is rather difficult to compute the integration in the relation (8) analytically due to the pair-distribution 
function g{f, r'). 

Since the HF approximation has already contained the contribution from the exchange effect between two electrons with 
the same spins, the correlation energy may be obtained by subtracting the HF pair-distribution function g^f(r,r') from the 
exchange-correlation energy E xc 28 ' 32 , 

2 aa > J \r-r'\ 
where p c aa , is the spin-dependent correlation hole and given by 

p% a , (f, ?) - Pa , {?)[g aa , (r, ?) - 5 aa ,g^ F , (f, ?)] . (12) 
where g^(r, r') is the pair-distribution function under the HF approximation, and is given by 

g%{?, ?) = f d\g™(r, ?■ A). (13) 
Jo 



The difference between the two pair-distribution functions (g — g HF ) has to do with the electron correlation. The expression 
(11) now can be regarded as an integral between the 7r electron at the position f and the charge cloud of the spin-dependent 
correlation hole around the position r' of the other n electron. Thus the expression (11) can be rewritten as 

E ° = \ E / d 3 rd 3 r' Pa (r)v(r- r')p c aa , (r, ?) (14) 

crcr' 

where v(r-r') = , e ->, , the Coulomb interaction between two electrons at the positions r and r'. 

\r-r'\ 

In order to complete the calculation of the integration of E c , an approximation is needed to simplify the expression for E c . 
In the study of polyacetylene (PA) oligomer, two-electron interaction integral was approximated by abstracting the Coulomb 
interaction v(f — r') out of the integrand of the many-centered Coulomb integral, and the results showed that the approximation 
is reasonable 33 . According to that spirit, we may take v(f — r') out of the integrand in (14). For simplicity, in this study 
only consider the nearest-neighbor electron-electron interaction between two adjacent carbon atom sites are considered, that is, 
v(r — r') = Xv, where A is due to the replacement of the charge e 2 in the Coulomb interaction. In this way the expression for 
E c may become 

^ aa' 
aa' 

= yE / d\d\'p a (r)p al (?)g aa ,{ry)-Y, f d 3 rd 3 r> ' p a {r)p a {?)g^ F \f,P)} (15) 

aa' J a J 

Remarkably, it is not easy to calculate it without further approximation. According to the sum rule (10), we have 

Y, I d 3 rd 3 r'p a (r)p a ,(P)g^(r,P) = 1 £ f d\d\' p a (r)p a ,(?) I d 3 rp a (^ (16) 

aa' J aa' J a J 

Then the expression for E c can be reduced to 
Xv, 1 



Xv 

T 



E c = 

2 



~ E / + \ E / d 3 rd 3 r' Pa (r)Pa>(P) J2 I d 3 rd 3 r' p a {r)p a {?)g^ F (r, ?)]. (17) 



Now the remaining central task is to calculate the Hartree-Fock pair-distribution function g HF {r 1 r l ). From (7), we have 

HF p = Ejfc(i#fc) < *HF\fllMk\*HF > \Mrl\ 2 \Mr')\ 2 

E£ < *HF\%i\*HF > \Mr)\ 2 J2l < QhfWQhf > \Mr')\ 2 ' 

where \&hf > is the HF ground state. 

Inserting (3) and (18) into the expression of E c , it yields 



E ° = t [_ ^E^ + ^ E «1-E<^f|«i^>]. (19) 



The third term in the above expression can be evaluated and it equals 

X) < ®HF\puPkk\$HF > = < ®HF\4 a a ia al lT aka\$HF > 

a aa' 

= E[< a ™ a ™ ><: a la a ka > - < 4a a ka >< a\a a ™ >\ 
a 

= E^fc-(PSb) 2 ]- (20) 

a 

Here the symmetry is used that < a\ a aka >—< o/ kcr ai a >. Finally the correlation energy obtained is 

E * - yI-^E^+E^) 2 ]- ( 21 ) 

ia ika 

This expression (21) is the central result in this paper and it gives the correlation energy for an open or closed one-dimensional 
7r electronic system in real space. The correlation energy per electron is e c = E c /N, N is the total numbers of the carbon atoms 
in a 7r electron system. For a half filled system, the number of atomic sites and the number of electrons are equal. This relation 
tells us that when the averages of the bond charge and the electron density at the site i are known, the correlation energy can 
be evaluated. 



4 Calculation and results 



In this study, the main consideration is about the correlation energy of the ir electron systems. For the ir conjugated polymers, 
the Hamiltonian of the system is the SSH-type Hamiltonian, H , plus the electron-electron interaction term, 

H = H ° + \ E ~ ^kW«,v a ^' ( 22 ) 



#o = ^ M^W + ^0 (23) 

ij(T(T f 

where a| cr (a :)(T ) is the creation (annihilation) operator of an ir electron at the site i (j) with spin a. v(r ~ r') is the electron- 
electron interaction, and r (r') means the position vector of an ir electron at the site i (j). tij is the hopping term. For 
an bond-alternated chain, tj,j+i = to + (— l) l 5<o with to being the hopping integral without dimerization and Sto being the 
magnitude of the dimerization due to Peierls transition. 

For the one-dimensional ir electron conjugated polymers with N carbon atoms such as the bond-alternated chain PA, in the 
bond order wave (BOW) phase, the average charge density at the site i p^ = 1/2, and the average of the bond charge density 
Pik = Pki =< a \a a k<j > ( ncrc k = i + Dropping the spin index because of pf k = p a lk = p lk , we have p ii+1 = p + {-l) l 5p. 
Then (21) becomes 

„ Xv r N n 9 

E c = yhy + 2E^ + i] 

i 

\ AT N / 2 

Xv . TV n \ . , 9 m ,„.s 

= Y^ _ T + 2 E^2m-l,2m + P2m,2m+l)]- (24) 
m— 1 

Since p2 m -i,2m = (P~dp) 2 = p 2 — 2p~5p + (5p) 2 and p\ m ^m+\ ~ (p + ^p) 2 = P 2 + 2pSp + (Sp) 2 , so the correlation energy for PA 
may be expressed as 

\v 1 

E c - -N-[--2(p 2 + {8p) 2 )]. (25) 

where N is the carbon atom numbers in the system. The correlation energy per 7r electron is then given by the integral over 
the parameter A from to 1, 

e c = - V -{j^d\-2j\\(p(\)) 2 +\(5p(\)f]d\}. (26) 
Here p(X) and 5p(X) are given by the first and second elliptic integrals 34 

KX) = n(l-\(\r) [E{y/l ~ Zm *W 2 K(y/l~zW 2 )] (27) 

SpW = ^^^ [^(v/1-^A) 2 ) - S(Vl-^A) 2 )] (28) 

where z(X) = St/t and the parameters t and St are determined by the relations 

t = t + Xvpo (29) 



St = Sto + XvSpo (30) 

where po and Spo are the average charge density and the change of the charge density without the Coulomb interactions. 

For the PA, t =2.5 eV and St = 2au = 0.269 eV, where a = 4.1 cV/A and the dimerization u = 0.03 281 27 . In 
calculation, the lang-range interaction v is 2.4 eV 27 . Table 1 lists the values for different A when using the the relations (29) 
and (30) to calculate the average of the bond charge density p and its change Sp from the relations (27) and (28). Putting 
the resulting p and Sp into the (26) and integrating over A from to 1, the correlation energy is then obtained. In Table 1 
I(X) = (p(X)) 2 + (Sp(X)) 2 , and I = 2 Jq[X(p(X)) 2 + X(Sp(X)) 2 ]dX. Table 2 lists the values of the correlation energies for PE 
and PA for contributions from the on-site Coulomb interaction and the nearest-neighbor Coulomb interaction, and also lists the 
band energy gap E' in this study and that from Ref.27. Fig.l shows the correlation energy vs the dimerization Sto- 

To see the influence of the parameter A on the the bond charge density and finally on the correlation energy, I calculate the 
p(X) and Sp(X) (see the lines from 2nd to 11th in the Table 1). It is seen from the Table 1 that the differences about the various 
quantities when A = and when A = 1 are very small: p(0) - p(l) = 0.0026, \5p(0) - Sp{l)\ = 0.0200, |7(0) - = 0.0024. 



Table 1: To see the influence of the parameter A on the the bond charge density and finally on the correlation energy, the p(X) 
and Sp(X) are calculated (see the lines from 2nd to 11th in the Table). It is seen that the differences about the various quantities 
when A = and when A = 1 are very small, to = 2.5, u = 0.0328, a = A.leV/ A and v = 2.4 for PA. The unit is eV. 
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In addition, the differences between 7 and I(X) are also small: |7 - 7(0) | = 0.0017 eV, |7 - 7(1) | = 0.0007 eV, and |e c - e c (A = 
0)| = 0.0020 eV, |e c - e c (A = 1)| = 0.0009 eV. 

The treatment of the pair-distribution function (9), g(f, r 1 ) = J Q dXg(f,r'; A) is as follows. In the formula (9), Coulomb 
interaction e 2 in the Coulomb interaction v(f — r') is replaced by Ae 2 , A changes from to 1. When v(f — r') appears and 
changes, the electron wave functions follows the changes. Then the electron density distribution p{r) changes, and then pij 
changes, which may be expressed as p(X) and Sp(X) (see (27) and (28)). This can be seen through the expressions (29) and 
(30). Therefore the calculation of g(f, r") is realized actually through the calculation of p(f)(X). Table 1 shows the results of 
the correlation energies with integration about the parameter A (see the first line in the Table 1). 

Note that the bond charge density pfj always is less than half and about 0.3 or so, therefore the second term is less than the 
first term in the formula (21) and the correlation energy is negative. The electron systems to which the formula is suitable is 
supposed to be 7r conjugated polymers with long chain (N is very large) with half filled band in the ground state. For a small 
molecule system such as 77 2 , etc, and those without ir electrons, the formula is not suitable because the molecule like 77 2 is 
covalent molecules where the charge density gathers between two atoms. Another reason is that the approximation (15) may 
bring a bigger error about the electron-electron interaction integral if it is applied to small molecules like Hi etc. It is also not 
suitable to the hydrocarbon such as methane molecule that has no 7r electrons. In derivation of the formula, there is no the 
excited states to be dealt with. 
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Figure 1: e c changes with increasing dimerization 5to (= 2auo) under v = 2.4 eV according to the formula (26). Here v in e c (v) 
means the nearest-neighbor Coulomb interaction. 



Table 2: Correlation energies e c from mtij and n^n^ for PA and PE. Here e c (V) means contribution from the nearest-neighbor 
Coulomb interaction, e c (U) means contribution from the on-site Hubbard interaction, and e c (U+V)means contribution from 
both the on-site Hubbard interaction and the nearest-neighbor Coulomb interaction. E' g (= E g + e c ) and Eg are the band energy 
gap. The unit is cV. 

Systems PE PA 



E'g 1.7832* 

E' g ' 1.8 d 

e c (V) -0.1725 a -0.1696* 

e c (V) -0.1567 6 

e c (V) -0.0100 c 

e c (U) -0.5760 a 

e c (U) -0.7885 6 

e c (U) -0.6208 c 

e c (U+V) -0.75 a 

e c (U+V) -0.6319 6 

e c (U+V) -0.6308 c 



a refers to the reference [22] 
b refers to the reference [23] 
c refers to the reference [24] 
d refers to the reference [27] 
* refers to the present study 



5 Discussion 

In the process of derivation for the formula, there are two approximations to be used. One is the tight-bind approximation where 
the pair-distribution function g(f, r') only involves diagonal elements pa of charge density. Other is that the electron-electron 
interaction v(r — r') is drawn out of the integrand of the integral expression for E c (See (15)), which was proved to be available in 
calculating the excitation of the conjugated polymer oligomer 33 . In this way, the correlation energy expression obtained is directly 
related to the diagonal site charge density p^ —< a^a^ > and square of the bond charge density p a ik =< a} ia aj a > (i ^ k) that 
can be calculated directly from the elliptic integrals (27) and (28). This is different from the variational methods 22 ' 23 ' 24 , and 
also from the the density matrix renormalization group (DMRG) 34 . 

If the overlap effect of the 7r electrons at adjacent sites is included, the tight-bind approximation is invalid and the pair- 
distribution function g(f,r') will have both the diagonal terms pa and the off-diagonal terms pij of the charge density. In this 
case, the electron charge density pa on the sites will be less than half and the bond charge density pij (i ^ j) will be enhanced 
a little. It may be predictable that when the overlap effects of the it electron waves are considered, the correlation energies 
will become smaller. Moreover, when the overlap integrals are considered, we may discuss the contributions of the correlation 
energy from the off-diagonal electron-electron interactions 35 . 

It is seen from Table 2 that for an infinite polymer PE, the correlation energy was obtained and was e c (U + V) = —0.75 
cV 22 . According to Ref.22, 77% the correlation energy comes from the contribution of the operator n^n^. Thus the remaining 
23% of the correlation energy comes from the contribution of the long-range Coulomb interaction Vij (i ^ j) or equivalently 
from the contribution of the operator n,nj (i ^ j). From this, we may have e c (V) — 0.23 x (—0.75) = —0.1725 eV. The present 
correlation energy (e c (V)) comes from the contribution of the nearest-neighbor Coulomb interaction (v = 2.4 eV) between two w 
electrons and is —0.1696 eV. In the Ref.23, the long-range interaction V (here V = V^+i = v) was included within an "effective 
U" after using some approximation, and the correlation energy e c (U + V) both including the on-site Hubbard interaction U and 
the long-range Coulomb interaction V and e c (U) were obtained for PA, then the correlation energy e c (V) from the contribution 
of the long-range Coulomb interaction V may be obtained from the difference between e c (U + V) and e c {U) and it yielded 
—0.1567 eV with t — 2.9 eV. This value is smaller than the result of the present study. The Ref.24 used the local approach 36 
to calculate the correlation energies. However, the difference between e c (U + V) and e c (U) was very small and e c (V) = —0.0100 
cV with t = 2.5 eV in Ref.24, and this result seems so small. 

It is also seen from Table 2 that the present result of the correlation energy e c (V) for PA is little smaller than that for PE. 
Both PE and PA includes sp 3 hybridization. In the calculation of the correlation energy for PE 22 , except it electrons between 
two nearest-neighbor carbon atoms , factors from the different w bonds were also taken in numerical computation, so the more 
correlation effects were included in e c (U + V) for PE. In addition, the present result e c (V) = —0.1696 eV is little larger compared 
with e c (V) = —0.1567 eV from Ref.23. This may be caused by the approximation where the overlap effect between two adjacent 
7r electronic wave functions is omitted. When the overlap effect is considered in calculation, the bond charge density will be 
larger and the site charge density will be less than half, then the result will become smaller. 

It is seen from (26) that the correlation energy is an even function of dimerization. The trend of the curve in Fig.l is kind 
of quadratic but not complete because there are also the dimerization parameter z in the denominators in the formula. This 
point can be seen by the following way. Because z <C 1, if we replace z in the denominator in the formula by zq(= St /t n ), then 



the correlation energy e c is approximately proportional to — 0.125v + ^ + Jpr-z 2 , where A and B are the integral constants. 
At present, there is no similar curve to compare. We may compare the curve in Fig.l with those in Refs.24 and 27. Because 
the correlation energy (26) do not contain the contribution form the on-site Hubbard interaction U, the curve in Fig.l is not 
completely like those in Fig. 2 in Ref.24 and those in Fig.l in Ref.27 where both U and the long-range interaction were included. 
Nevertheless, it could be found here that the curve in Fig.l still have some similar trend with them when the dimerization value 
is bigger. Though the curve in Fig. 2 in Ref.27 was about the ground state energy vs the dimerization, the curve also reflects 
some information about the correlation energy vs the dimerization because the ground state energy contains the correlation 
energy contribution in Ref.27. 

It is all known that when screening is weak or normal in the it electronic conjugated polymers, the electron-electron interaction 
increases the dimerization and band energy gap 35 . That is, v increases Sto and E g . In equilibrium state and rigid background, 
the band energy gap E g of PA is given by ASt = A(8to + vSpo). Under the electron-electron interaction v, the average bond charge 
density po decreases slightly with increasing Sto and the fluctuation of the bond charge density Spo increases with increasing 
Sto. The decrease of po causes the bandwidth (see t = to + vp ) to diminish, and the increase of 5po make E g increase and 
at the same time cause e c (see (25)) to decrease. Then it can be seen from these that v and 5po are two opposite factors to 
the correlation energy: v is in favor of the correlation energy but Spo is a disadvantage to the correlation energy in this study. 
That is to say, on the one hand, v makes Sto an d accordingly dpo increase and then the band energy gap E g increase, on the 
other hand, the fluctuation 8po will cause the correlation energy e c to decrease. As a result, these two opposite influence makes 
E g decrease from E g = ASt to E' = E g + e c . When Sto = 0.269 eV and the corresponding bond charge density Spo = 0.0903, 
we have E g = 1.9429 eV and e c = —0.1696 eV. Therefore the band energy gap E' containing the nearest-neighbor Coulomb 
interaction correlation effect becomes 1.7832 eV. This value of the band energy gap is close to E g (1.8 eV) obtained by Ab initio 
computation by author in Ref.27 where the screened interaction was used. These are in qualitative agreement with experiment. 
In addition, E' g increase with increasing dimerization because the correlation energy e c decreases with dimerization, which is 
consistent with relation of the band gap and the dimerization. 

To my knowledge, there is no similar expressions for the correlation energy only from the nearest-neighbor Coulomb inter- 
action that exists in a simple form at present. In the next study, the more electron-electron interaction terms in the long-range 
Coulomb interaction will be considered. In addition, although the DMRG is a strong tool to deal with the correlation problems 
in a many-particle electron system, it is basically a complicated numerical calculation method but not is an analytical expression. 

In summary, under the approximation (15) and the tight-bind approximation, a formula (21) of the correlation energy 
for the long-range (nearest-neighbor) Coulomb interaction v for the conjugated polymers is obtained with the rigid backbone 
background. Although it is simple, it is direct and effective and easily operational in comparison with other highly involved 
numerical computation methods including DMRG. The computational result for the correlation energy for PA is available 
compared to those for PA and PE in different methods 22 ' 23 . The band energy gap E' g containing the correlation effect is close 
to that by ab initio method containing the screening interaction 27 . Because there are no constraints to the systems in the 
process of the deduction, this relation may be applied to the various n electron systems such as Ceo, benzene rings, and carbon 
nanotubes, etc. 
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